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DEFINITION OF SYMBOLS 


Symbol 

( 0 ) 


4 > 


n 


H 


(o) 

s (o) 

J n 


n 

H 

X 


(i) 


n 


H 


(i) 


n 


n 

(i) 

n 

(i) 

n' n 

s ( :> 

n' n 


E 


H 


Definition 

Unperturbed wave function for the system in the n state 

Unperturbed Hamiltonian 
Unperturbed energy eigenvalue 

Wave function for the perturbed system 

Hamiltonian for the perturbed system 

Expansion parameter, such that A < 1 

For i + 0, these are the wave function corrections for 
the perturbed system 

For i + 0, these are the perturbation terms in the 
Hamiltonian 

Energy for the perturbed system 

For i + 0, these are the energy corrections for the 
perturbed system 

For i + 0, these are the matrix elements obtained using 
the unperturbed wave functions 

For i # 0, these are the matrix elements defining the 
unitary matrix 


iv 



APPLICATION OF UN IT ARY TRANSFORMATIONS 
TO FOURTH ORDER PERTURBATION THEORY 


SUMMARY 


Perturbation theory is used to obtain energy corrections to the energy 
of a nondegenerate stationary system. The unitary transformation method is 
used to obtain the energy corrections to fourth order. In the first section, a 
short introduction is given to the usual perturbation techniques. In subsequent 
sections, the higher order corrections are given in a somewhat detailed manner. 
The purpose of this approach is to point out the difficulties encountered in 
using a formalism to obtain useful formulae suitable for numerical calculations. 


I. INTRODUCTION 


Quantum mechanical problems which can be solved exactly take on the 
simple form 


„<»> 





( 1 ) 


In this discussion the nondegenerate, stationary type of problem is assumed. 
The above equation then applies to the unperturbed system, where H^°) is the 
Hamiltonian. If small perturbations are introduced into the system, the 
Hamiltonian is usually made up of terms in addition to H^°). In this case 
equation (1) can no longer be written in closed form, and 'k is usually found by 
means of a series expansion made up of successive approximations which are 
denoted by 0 , thus 


<f? = + A 

n r n r n 


+ X 2 + 
r n 


A < 1 


( 2 ) 


Here A is a small parameter which goes to zero when the perturbation is re- 
moved. The above wave function then corresponds to a modified Hamiltonian 

+ A 2 H (2) + ... 


H = H (0) + A H (1) 


( 3 ) 



In the above expression, the Hamiltonian terms due to the perturbation are 
separated into various factors depending on their magnitude, the larger terms 
appearing first. In the above formalism, tp^ is the exact solution to the 
wave equation when only appears in the Hamiltonian. If H is separated 
as indicated, the energy eigenvalue may also be expressed as a series expan- 
sion as follows: 

E = E (0) + X E (1) + X 2 E {2) + ... . (4) 

n n n n v ; 

If one substitutes the above values in Sehroedinger's equation, one obtains, 


(H (0) + A H (1) + ...) (^ (0) + A ^ 


+ ...) = (E (0) + X E (i) + ...) 
' n n f 


X (/ 0) + X ip^^ + ...). 

' Y n n ' 

At this point one needs to decide how far to carry on the approximation. 

A first order approximation includes terms involving X to the first 
power, all other terms are neglected. Thus to first order, 

(H (0) + X H (1) ) (/ 0) + X 4> (1) ) = (E (0) + X E (1) ) (j p (0) + X ip (i) ) 
n n n n ' Y n Y n ' 


or 


/"> -E ( "> /”>) (H (l > ♦< 0 > + H (,) *<‘> - E ( ‘> /«> 

n n r n ' Y n Y n n Y n 


- E^ ip ^ ) = 0 
n n ’ 


Using equation (1) , the first term vanishes. In order for the rest to be true, 
the X coefficient must vanish. From the last term above, one gets 


(H ( °* - E (0) ) + (H^ - E (l) ) = 0. 

n ' r n n ’ r n 

Operating on the above equation as follows , one obtains 

(^ 0) | H (0) - E (0) lip (l) ) + (ip (0) IH (1) - E (l) l ip (0) ) = 0 , 
m n n ' vv m n r n ' 


2 



(E 


(o) „(0). 


m 


- E w ) 8 
n ' mn 


( ^ 0) |H (1) |^ 0) ) = E (1) 6 


r n 


n 


mn 


Since one deals with orthonormal functions, the above simplifies to 

E (1) = (/ 0) |H (1) | i =H (1) . 
n r n n nn 

In integral form one gets the usual form 

E (1) = f i H (l) / 0) dT. 
n J r n T n 


( 5 ) 


In this manner one may obtain all the terms in equation (4) which correspond 
to the corrections to the energy E n of the system. For higher order corrections, 
the algebra rapidly becomes too complex. 
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II. APPLICATION OF UNITARY TRANS FORMAT IONS 


For higher approximations an alternate method [ 1] is much easier to 
develop. One simply connects with ^ through a unitary transformation and 

uses a representation in which forms a diagonal matrix such that 




= = E (0) 6 

n' n n n' n 


With the above requirements, a unitary matrix may be defined such that its 
elements are given by 


U 


n' n 



11* ) 

n 


( 6 ) 


3 



Using such a matrix, the Hamiltonian of the system in question may be diago- 
nalized by taking 


H' = U H U" 1 . (?) 

After H is diagonalized, the matrix elements of H' are just the required energy 
corrections. 

In order to see how the matrix U is developed, one must recall that the 
tp^'s are the basis states corresponding to the unperturbed system, and that 
the new 4>'s form a new basis corresponding to H. Thus, if one lets 


* 

n 


N 

S u* 

3=1 nJ 



then one can also write 


(<k , | ) = 5 

n' n n' n 


( 8 ) 


(9) 


Substituting equation (8) in equation (9), one obtains 
N N 


y y u u*. 

i-i n'i m 

j=i i=l 


. ( 1 ^°^) =5 
ni j l n’n 


Simplifying by using the orthonormality condition of the ip^ 's, one obtains 


N 


y u u*. = 6 , 

t-j n'i ni n'n 

j=l 

or by using the <5 condition explicitly 


N 


y u . u*. = i. 

U nj nj 


j=l 


Writing the above in matrix form, one gets 
U U t = 1, 


which implies that 


4 



u + -ir‘ 


This condition says that U has an inverse which is equal to its Hermitian 
conjugate. This U, in turn is a unitary matrix capable of diagonalizing the 
Hamiltonian H which amounts to solving the problem. 

If this method is used in obtaining the energy corrections, the unitary 
matrix may be written in exponential form. Thus, one can define U as 

U = e iS , i 2 = -l. (10) 

S must be a Hermitian matrix (S = S^~) in this case. With the above definition 
for U one can easily see that S may be expanded in a series of the form 

s = x s (1) +a. 2 s (2) + ... . (11) 

This expansion is such that when X goes to zero, U becomes the unit matrix as 
necessary. Before proceeding it is desirable to expand the newly diagonalized 
Hamiltonian in a similar form, thus 


H' = E (0) + X E (1) + X 2 E (2) + . . . . (12) 

/ \ 

H' is expanded in terms of E' * 's because they now correspond to the desired 
energy corrections in the final result of equation (4) 

To obtain the diagonal elements of H';one starts by substituting equation 
(10) inequation (7). In addition, an arbitrary parameter (x) is introduced which 
will allow the eventual expansion of H 1 in a Taylor series. Hence, from equation 
(7) one obtains 


H’ =H(x) = e lxS H e lxS 


(13) 


Expanding H (x) , one gets 


/ dH(x) 1 

1 

d 2 H(j£l| 

1 d * J 

X + 2 ? 
x=o 2 - 

l d X 2 1 


X^ + 
x = 0 


When x goes to zero, H’ is still well defined, that is 
H’ = H(0) = H. 


(14) 


(15) 
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In order to find the terms of H(x) , one uses equation (13) to get the indicated 
derivatives. Remembering that S is independent of x, one gets 


- = 1 S e lxS H e-‘* S + e‘* S H e' 1 * 8 (-IS) 

= i(SH(x) - H(x) S) *i[S, H(x) ] . 

In the last line the bracket notation is introduced to make things simplier later 
on. Using the above notation the second derivative is obtained as follows: 


d 2 H(x) _ d 

dx 2 dx 


(i [S, H(x) ] ) - i [S, 


d. H(x ) 

dx 


= i 2 [S, [S, H(x) ] ] . 


In general [2] , the a derivative is given by 




dx 


= [S, [S, ... a , [S, H(x)] ... a ,]). 


Equation (14) becomes 


H( X ) =H(0) +i (S, H(0) ] x + £T[S, [S, H(0) ] ] X 2 + . . . • 

When one lets x = 1 (its value is arbitrary) and uses equation (15) the above 
equation becomes 

H' = H(1) =H + i[S, H] + “7” [ S, [S, H] J + ... . 

By means of equation (13) , one sees that H(i) is just H\ To fourth order, H' 
is given by 


H’ = H + i[S, H] +— [S, [S, H] ] +— [S, [S, [S, H] ] ] 


+ -£7 [S, [,S, [S, [S, H)]]]. 


(16) 
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In subsequent discussions higher order terms will be omitted. Using equations 
(3) and (11), one obtains for H* the following expression: 



H' = (H (0) + X H (1) + X 2 H (2) + . . . ) 

+ i [ (X S (1) + X 2 S (2) +...), (H (0) + X H (1) + X 2 H (2) +...)] 

+ 7 T[(XS (1) + X 2 S (2) + ...), [(\S (1) +A 2 S (2) + ...), 

£Jw 

(H (0) + X H (l) + X 2 H* Z) +...)]] 

+ ^7tS, fS, [S, H]]] + ^7 [S, [S, [S, [S, H] ] ] ] . t 16 - 1 ) 

Expanding the first bracket in H* one gets 

{[(X S (l) + X 2 S (2) + X 3 S (3 > +X 4 S (4) ), H (0) ] 

+ [ (X S (1 > +X 2 S (2) +X 3 S (3) + X 4 S (4) ), XH (1) ] 

+ t(x s (1) + x 2 s (2) + x 3 s< 3) + X 4 S (4) ), X 2 H (2) ] 

+ [(X s (1 > +x 2 s (2) +x 3 s (3) +X 4 s (4) ), X 3 H (3) ] 

+ [ (X s (1) + X 2 s (2) + X 3 S (3) + X 4 S (4) ), X 4 H (4) ]} . 

By neglecting all terms in which X is raised to a power greater than four, one 
simplifies the above expression and obtains 

{[XS (1 \ H (0) ] +[X 2 S (2) , H (0) ] +[X 3 S (3) , H (0) ] + [X 4 S (4) , H (0) ] 

+ [X S (1) , X H (l) ] + [X 2 S (2) , X H (1) ] + [X 3 S (3) , X H (1) ] 

+ [X S (1) , X 2 H (2) ] + [X 2 S (2) , X 2 H* 2) ] (Bi) 

+ [X S (l) , X 3 H (3) ]> 

the above expression is written such that it shows how each term is obtained. 

The first bracket has been evaluated in detail because it is used to evaluate sub- 
sequent terms. This can be seen by reference to equation ( 16) . Thus in the 
calculation of the next bracket in equation (16) , use is made of the results obtained 



for the first bracket as they appear in expression (Bl) . Expanding, one gets for 

the next term of equation (16. 1) , [ (X + X 2 + ...),{[ S, H] }] , the 
following: 

{[X S {1 \ {[XS (1) , H (0) ]+ [X 2 S (2) , H (0) ]+ [X 3 S (3) , H (0) ] + [X 4 S (4) , H (0) ]}] 

+■ [X 2 {[X ] + [X 2 S^, H^] + [X 3 S^ 3 *, H^] + [X 4 , H^]}j 

f [X 3 s (3) , {[XS (l) , H (0) ] +[X 2 S (2) , H (0) ] +[X 3 S (3) , H (0) ] + [X 4 S (4) , H (0) ]}} 

+ [X 4 S (4 \ {[XS (1) , H (0) ] +fX 2 S (2) , H (0) ] + [X 3 S (3) , H (o) 3 + [X 4 S (4) , H (o) ]}] 

+ [X S (1) , {[X S (l \ X H (l) ] + [X 2 S (2) , X H (l) ] + [X 3 S (3) , X H (l) ]}] 

+ [X 2 S (2 \ {[X S (l) , X H U) ] + [X 2 S (2) , X H (l) ] + [X 3 S (3) , X H (l) ]}] 

+ [X 3 S (3) , {[X S (l) , X H (1) ] + [X 2 S ( \ X H (1) ] + [X 3 S (3) , X H (1) ]}] 

+ [X 4 S (4) , {[X S (1) , X H (1) ] + [X 2 S (2) , X H (1) ] + [X 3 S (3) , X H (l) ]}] 

+ [X S (l) , {[X S (l) , X 2 H (2) ] + [X 2 S (2) , X 2 H (2) ]}] 

+ [X 2 S (2) , {[X s (l) , X 2 H (2) ] + [X 2 S (2) , X 2 H (2) ]>] 

-l- [X 3 S (3 \ {[X s (1) , X 2 H (2) ] + [X 2 S (2) , X 2 H (2) ]}] 

+ [X 4 S (4 >, {[X S (1) , X 2 H (2) ] + [X 2 s (2) , X 2 H (2) ]}j 

+ [X S (1 \ [X S U) , X 3 H (3) ]] 

+ [X 2 S^, [X , X 3 H^] ] 

+ [X 3 s (3) , [X s (l) , X 3 H (3) ]] 

+ [X 4 S (4) , [X s (l \ X 3 H (s) ]]} . 
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By neglecting terms having \ n , n > 4; one gets for the second bracket in equation 
( 16. 1) the following expression: 


{[X S (l) , [A S (l) : 


H (0) ] ] 


n 2 tj(°) 


(l) 


+ [A S'*', [ A d S'", H' ; ] ] + [A S'", [A d S'", H' ' ] 


,3 H (3) tj(°) ■ 


+ [A 2 S< 2 ), [AS< 1 >, H (o) ]] + [ A 2 S^, [A 2 S^, H ^ °^ ] ] 

+ [A 3 S (3) , [A S (1) , H (0) ]] 

+ [A S (1) , [ A S (1) , AH (1) ]] + [A S (1) , [A 2 S (2) , A H (1) ]] 


(B2) 


+ [A 2 S (2) , [A S (1) , A H (1) ]] 


+ [A S (l) , [ A S (1) , A 2 H (2) ]]} . 


Clearly the last factors on each one of the lines in expression (Bl) did not con- 
tribute to expression (B2) . Hence, in the calculation of the third bracket in 
equation (16) , the last factors on each one of the lines of expression (B2) will 
be left out. With this in mind, the third bracket of equation (16) is obtained as 
follows: 

By taking 

f(AS (1) +A 2 S (2) +...), { [S, IS, H ] ]} ] , 
and by substituting the appropriate terms in expression (B2) , one obtains 


{[AS 

(l) 

{[X S (l) , 

[A S (l) , 

H (o) ]] + 

[A S 

(1) 

> 

[A 2 S 

<2> , H 

(0) ]]}] 

+ [A 2 

s< 2 > 

, {[A S U) 

,[xbW 

, H (0) ] ] 

+ [ A 

g (i) 

. [x 2 

s (2> , 

H (0) ]]}] 

+ [A 3 

S<3> 

, {[A S (l) 

, [AS (1) 

, H (0) ] ] 

+ [A 

s (i) 

. [X 2 

s (!> , 

H (0) ]]}] 

+ [A 4 

g(l) 

, { [A S (l) 

, [ A S (l) 

, H (0) ] ] 

+ [X 

s< 1 ) 

, [X 2 

s (2) . 

H (0) ] ]} ] 

+ [A : 

3 ( *>, 

[A 2 S (2) , 

[X s (1) . 

H (0) ] ]] 







+ [A 2 S (2) , [A 2 S (2 >, [X s (l) , H (o) ]]] 

+ [A 3 S (3) , [ A 2 S (2) , [A S (l) , H (0) ]]] 

+ [A 4 S (4 \ [A 2 S (2) , [A S (l) , H (0) ]]] 

+ [A S (1) , [A S (l \ [A S (l) , A H (l) ]]] 

+ [A 2 S (z) , [A S (l) , [A S (l) , A H (l) ]]] 

+ [A 3 S (3) , [A S (l) , [A S (l) , A H (1) ]]] 

+ [A 4 S (4) , [A S (l) , [A S (1) , AH (l) ]]]} . 

Several terms could have been omitted from the above expression but were re- 
tained only for symmetry. Excluding the terms having A n for n > 4, one derives 
for the third bracket of equation ( 16) the following expression: 

{[AS (1) , [AS (1) , [AS (1) , +[AS (1 \ [AS (1) , [A 2 S (2 \ H (0) ]]] 

+ [A 2 S (2) , [AS (1) , [AS (I) , 

+ [AS (1) , [A 2 S (2) , [AS^\ H (0) ]]] 


+ [A S (1) , [A s (1) , [A S (1) , A . 


(B3) 


For reasons previously stated, only the first term in expression (B3) contri- 
butes to the fourth bracket. From equation (16) and expression (B3) there 
results the following expression for the last bracket in equation ( 16) : 

[(A S (l) + A 2 S (z) +...), [A S (l) , [A S (1) , [A S (l) , H (0) ]]]]. 

If this term is expanded, only the first factor is retained, i. e. , 


[AS (l) , [AS (1) , [AS (i) , [AS (1) , H (0) ]]]] 


(B4) 
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Incorporating the expressions found (Bi, B2 f B3, B4) for the various brackets 
in equation (16.1) H’.becomes 


H’ = + A H* 1 * + A 2 H* 2 * + A 3 H* 3 * + A 4 

+ i{A (S (l) , H (0) ] + A 2 ([S (2) , H (0) ] + [S (l) , H (l) ]) 

+ A 3 ([S (3) , H (0) ] +[S (2) , H (l) ] +[S (1) , H (2) ]) 

+ A 4 ([S (4) , H (0) ] + ( S <3 \ H (l) ] + [S (2) , H (2) ] + [S (l \ H (3) ] )} 

+ ~ (a 2 [S (1) , [S (l) . h (0) ]] + a 3 ([S (l) , [S (2) , h (0) ]] 

[S (1) , H (o) ]] +[S (1) , [S (l) , H (l) ]]) 

+• A 4 ([S (l) , [S (3) , H (0) )l +[S (2) , [S (2 \ H (0) ]] +[S <3) , [S (l) , H (o) ]] 
+ [S (1) , [S (2) , H (l) ]] +[S (2) , [S (l) , H (l) ]) +[S (1) , [S (l) , H (2) ]])} 

+ ^-{a 3 [s (1) , [s (l) , [S (1) . h ( 0 ) ]]]+a 4 (rs (1) , [s (l) , [s (2) , h (0) ]]] 

+ [S^, [S* 1 *, [S (l> , H* 0) ]]] + [S (l \[S* 2 \ rs (l \ H (o) ]]] 

+ [S (l) , [S (l) . [S (l) , H (l) ]]])} 

+ ~ o 4 r s(1) » [S (1) . [S (1) > I s(1) > h(0) ]]]]>- 


To solve for the EC*) elements one compares the equation above with 

equation (12) and simply equates coefficients of the various powers in A. Hence, 
one obtains the following equations: 


11 



(°) 

= H (0 > 

9 




(17) 

(i) 

= H<‘» 

+ i 

[s (1 >, 



(18) 

(2) 

= h (!) 

+ i 

([S (2) 

, «<“>[ 

+ [ s (l) > hU) ]} +^i[ S (l) , [S (l) , 

H (0) ] ]» 







(19) 

(3) 


+ i 

([S (3) 

, H (0) ] 

+ [S (2) , H (l) ] + [S (l) , H (a) ] ) 



[S W , H (0) n +[S'*', [S' 1 ', H' v '[ ] + [S'",[S'-',H'-'] I) 


,(l) 

( 20 ) 


,(2) rfi (l) „(0) 

+ ^-[s (l) , [S (l) , [S (l) , H (0) ]]], 

E (4) =H (4) + i ([S (4) , H (0) ] +[S (3) , H (1) ] +[S (2) , H (2) ] +[S (1) , H (3) l) 

+ 4([S (1) , [S (3) , H (0) ]] + [ S (2) , [S (2) , H (<)) ] ] +[S (3) , [S (l) , H (0) ] ] 

+ [S (1) , [S (2) , H (l) ]] +[S (2) , [ S< A > , H (l) ] ] + [S (1) ,[S (1) ,H (2) ]]) 

+ ^T ([S (l) , [S (i) , [S (2) , H (0) ] ] ] +[S (1) , [S (2) , [S (l) , H (0) ] ] ] 

o. 

+ [S (1) , [S (l) , [S (l) , H (l) ]]] +[S (2) , [S (1) ,[S (1) , H (0) ] ] ] ) 


+ ^[S (l) , [S (l) , [ S (l) , [S (l) , H (0) ]]]] . 


( 21 ) 


1 1 1. FIRST AND SECOND ORDER FORMULAS 


In the preceding section, formulas were derived for the first four cor- 
rections to the energy of a perturbed system. These equations will now be ex- 
pressed in terms of specific matrix elements. This is the final step in obtaining 
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actual energy corrections in analytical form, because the wave functions which 
are solutions to the unperturbed problem are used in evaluating these matrix 
elements. 

Since one uses a matrix representation in which is a diagonal matrix, 
equation (17) is simply an identity, that is, 


E 


8 = 

n n'n n'n 


Writing equation (18) in matrix form, one obtains* 

E (l) 6 , = + i V (S^,, - H^,, ) 

n n'n n'n ^ ' n'n" n"n n'n" n"n 


= + i (E (o) - E ( ^) S<;> 

n'n n n' n'n 


( 22 ) 


Since only the diagonal terms are desired, one can simplify the above expres- 
sion to 


E (i) 


n 


= H^ . 
nn 


As expected, this result agrees with equation (5), 

Starting from equation (19) , one obtains the second order formula as 
follows: Take equation (18) and solve for [S^, Substitute this value 

inequation (19). With this modification, E^ becomes 


^ 6 n'n = H n'!i + 1 ([s(2) ’ ^ + 2 lS ^ ’ + 2 ^ ^ > n'n * (23) 


The diagonal terms are obtained as follows: Let n' = n, then 


* Use is made of the matrix multiplication rules; 

c, . = / . b, . a.. and d. = Y V a, . b.. c. 

kj Y ki ij 1m Y j li ij j 


13 



E (z) = H (z) + i V {(S (2) t - H (0) , S^) 

n nn H nn' n'n nn' n'n 7 
n' 


+ £ (S (l) , H ( f ) - H (l) , S ( J ) ) + £ (S (l) , E (l) 6 , -E ( l ) 6 , S ( J ) )}. 
2 nn 1 n f n nn f n'n' 2 nn 1 n n*n n f nn f n f n ,J 


By replacing by and by recalling that is a diagonal matrix, 

one can reduce the above equation to this form: 


E ( 2 ) = h ( 2) + IV (S d) H (j) _ H (i) s (i) j 

n nn 2 ^ nn* n'n nn' n'n 7 
n' 


(24) 


-( 2 ) 


To solve for E 7 in terms of known matrix elements, one needs to obtain an 
n 

expression for . This can be accomplished by using the nondiagonal terms 


of equation (22). For n' # n, one obtains* 
0 = + i (E (0) - E^) 

n * n ' n r» * 7 n«r» 


n n 1 ' n T n 


or 


iH( ; ) 

,(i) n'n 

n'n 


n' + n. 


(E° - E 0 ,) 
' n n' 


(25) 


Since only the nondegenerate case is being treated no problems are encountered 
with denominators of the type encountered in equation (25) and subsequent work. 

When is substituted in equation (24) , one gets 


... , H (l) , H (1) . \ 

e (2) = h ( 2) + jt y 1 

n nn 2 “ 


n'#n 


nn' n'n 


\i (E° - E° t ) i (E° - E °)/ 
' n n' 7 v n' n 7 


and by rearranging terms, 

H (l) rr(l) 

e (2 >=h ( 2 > + £ nn ' n ' n 


n 


nn lj oo 
n'*n (E - E ,) 
n n' 7 


( 26 ) 


,0) 


* In the discussion that follows, the parenthesis around the zero in E' 7 will 
be dropped to simplify the notation. 1 
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This result can easily be verified by extension of the method outlined in the 

introduction. To solve for , one is required to know S^ f \ This will 

always be the case. In the next section, for instance, to solve for E^ ' one 

needs to know * 
n*n 


IV. THIRD ORDER FORMULAS 


So far, the formulas derived are not too difficult, and any method is 
equally suitable. For the third order corrections, the present method appears 

to have certain advantages. To solve for E^ ' , one substitutes values for 

[S^ 2 \ H^] and[S^, H^] in terms of E^ and E^ . From equations (18) 
and (19), one obtains 

|S<‘>, H<*>1= ( E<1> -H (1> > . 


[S 


( 2 ) 


(o) _ (e (2) -h (2 1i i (i) 
H i 2 S 




E (l) ] . 


After substituting these values into equation (20) and rearranging factors, 
equation (20) becomes 


E < 3 ) =h* 3 * + i ([S (3) , H (0) ] [S (2) , H (l) ] +^[S (1) , H (2) ] 


+ |[S (1) , E (2) ] +|[S (2) , E (l) ] ) 


+ -J~U s(l \ [S (l) , E (1) ] ] - [S (l) , [S (i) , H (l) ] ] ) . 
>( 3 ) , 


(27) 


In order to solve for E ’ in terms of known quantities, it is necessary to ob- 
(2) 

tain Taking the nondiagonal terms of equation (23) one obtains 
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Using 


0 = H n*n + i ([s(2)> r(0)] + irs (1) , h(1) i + h s(i) > e(1) ])„,„ - 


ivn 


= E^ 6 
n'n n n'n 


H (0) = E° 
nn n 


one obtains 


0 = + i (E° - E 0 ,) (E (l) - E ( ^) 

n'n n n' n'n 2 n n' ’ n'n 


+ L y /g(l) H (l) ttO) g(i) V 

2 U 1 n'n" n"n n'n" n"n ; 
n" 


After substituting for and solving for one obtains (for n'^n) 


S 


i H ( ^ 
(2) _ n'n 


n'n o o 
(E - E ) 
n n T/ 


. (H ( I\ - H (l) ) 

1 n 1 ]! n T n ! nn 7 

+ 

2 0 0 2 
(E - E ,) 

n n T/ 


1 Z 




n"*n,n' (E° - E°.) 

n n* 


o o 

E - E # 
n n f 


+ 


o o 
E n' “ E n " 1 


(28) 


( 2 ) 

The above solution for S , determines only the nondiagonal terms. The diagonal 

n'n 

terms are arbitrary and are not important in the energy corrections. These 

terms affect only the phase relations of the wave function corrections ip The 

(2\ n 
( * ) in the last term in S , indicates that, when this term is multiplied out, the 

n'n 

factor corresponding to n" = n is left out in the first term and the factor cor- 
responding to n" = n' is left out in the second term. Having obtained S^) , one 

( 3 \ 11 11 

proceeds to calculate E . One obtains the diagonal terms from equation (27). 

n / ■ \ 

By using the diagonal properties of E v , only the following terms remain 
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.(3) _ TT (3) 1 


E'” =H'-' +-i V f(S (!> , H ( ‘> -H (l) , S ( |>) + <s"' h7 - H'"’ s';')} 
n nn 2 Li lv nn' n'n nn' n'n 7 nn' n'n nn' n'n 7J 
n' 


,(l) „(*) „(2) o(i) 


-^([S (l) , [S (l) , H (l) ]] - [S (l) , [S (l) , E (1) ]]) 


nn 


In order to calculate the above, one must substitute for and 
term in equation (29) is obtained separately, as follows: 

The second term is 


n**n 


U | 

1 h (2) , h ( :* 

[ nn' n ! n 

h(1) , 

nn’ 

h (2) \ 

n'n j 

l 2 1 

1 (E° -E°) 
n T n 

( E° - 
n 

-E»,> 1 


(29) 

Each 


( H (l) , (H (1) - H ( )\) H ( )* 

1 nn T nn n f rr n’n 


H (1 * (H (1 * - H (l *) 

nn' n'n n'n' nn 


0 \2 


( E u - E u ) 
n' n 


0 \ 2 


( E u - E u ) 
n n' 


*2 

+ E E 1 

n'^n n"#n,n' 


r (1) (i) (i) 

) nn" n"n' n'n 

( 1 

1 \ 

l ( E° - E° ) 

v n' n 

U®-E® E«-E®/ 

n n M n T n M 

H ( 1 ) H (1) H (1) 

nn' n'n" n"n | 

1 1 

. . v 


(E° -E 0 ,) 
n n' 


Simplifying this term, one gets 


n' - n" 


E u - E u 
n n r 


i §{ 


m( 2 ) tt( 1 ) + Xj(i) H (l) (H (1) - H^) 

< H nn' n'n H nn' n'n* H nn' n'n < n'n' n nn > 


(E° - E°,) 
' n n' 7 


0 0 9 

(E - E ,) 2 
' n n' 7 
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(29.2) 


+ i Z Z 

n'^n n"^n, n' 


, H 0) H (l ) H (l > +H (i) H (i) H (l) \ 

* ' nn" n"n' n'n nn' n'n" n"n 7 


o o 
(E - E ) 
' n n' 


1 1 
+ 


l E° - E° E°,-E ° I 

\ n n ». n' n "/ 


The third term is 


i Z 


n'^n 


H (l) H (2) H (2) H (l) 

*' nn' n'n 


nn' n'n 


ii(E° - E°J i(E - E ) 
1 ' n n' n' n' 


-i Z 


(H (2) , H ( f ) + H (l * H (2) ) 
nn' n'n nn' n'n' 


n'#n 


(E° - E°,) 
n n' 


(29. 3) 


The fourth term is 


- -L S <‘> S (I) H (I) + ±S (l) H ( '> S ( ‘> -4- H<‘> S<‘> S<‘> 

6 


12 


12 


H • S w S' 


nn 


,( 1 ) 


Substituting for , one obtains 


Ti Z Z 


H (l) H d) H 0) 

nn T n T n n n M n 


LJ L—J o 0 ft o M ..V / vU 77 U \ / 77 U 77 U \ 

n'*n n"*n’ (E n - E^,) (E n , - E n „) n'/n n"*n <E„ - E^,) <E n - E n „) 


C° - E° ■ 6 

n T n 


I z z 


(1) (1) (1) 
H 11 . 11 

nn 1 n T n n n n n 


?0 _ 17O 


lz z 


h<‘>, h<;,> , h <!> 

nn" n"n' n'n 


12 U L-l 0 0 0 0 

n'*n n"*n' (E - E ,) (E - E „) 
' n n T/ ' n' n" 7 


The last factor was obtained by interchanging n' and n". Further simplification 
can be accomplished by breaking up the "1/6" term into two "1/ 12" terms and 
interchanging n" and n' in one of these terms before combining with the other two 
terms. By the above manipulation, the fourth term becomes 
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Ti l Z* 

n'^n n'^n, n' 


(H (1) , H ( }\, H ( ^ +H (1) „ H^) 

nn' n'n" n"n mi" n"n' n'n 


(E° - E° ) 
n n' 


o 

E - E , 
n n' 


E° - E° 
n' n' 


(29.4) 


The last term is 


_L a (0 a d) e (D _L a (D TC (i) a d) * ± *(*> s (i) s (i) 


112 


E' 


12 


E’ 


nn 


Substituting for one gets 


z 

n'=£n 


H (l) , E (l) 
1 nn 1 n’n n 

12 0 o o 

(E - E ,) 2 
' n n' 7 


, H ( ‘> H<;> E<:> , H (1 » H<:> E<‘) 

_1 nn* n'n n' _1_ nn* n’n n 

6 .„0 „0 .2 12 ,„0 _ „0 2 


(E - E .)' 
n n' 7 


(E - E ,)‘ 
n n' 7 


which, when combined, becomes 


H (l) , H ( J* (H ( J\-H (1) ) 
1^ y nn' n'n n'n' nn 

fi h 


n'^n 


( E - E ,)‘ 
n n' 


(29. 5) 


( 3) 

When all the above terms are collected, one obtains for E^ : 


e (s >=h ,3 >+ £ 

n T ^n 


n 


nn 


f(H (2) , H ( J ) + H (l) , H (2) ) H (l) , (H^, - H (1) ) 

1 nn' n'n nn' n'n 7 1 nn* n'n n'n' nn 

3 (E° - E°,) 2 

n n' 7 


o o 

(E - E ,) 
' n n' 7 


1 Z. Z 


❖ 


(l) (l) „(l) +h (D „(l) tt( i) v 

^nn' ^n'n" n"n nn" n"n' n'n 


n'^n n'^n, n' 


(E° - E° ) 
n n' 


o o 

E - E , 
n n' 


o „o 
E , - E ,, 
n' n' 


(30) 
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The preceding result may be simplified by modifying the last term so that the 
double sum will be symmetric in n ! as well as n". This is accomplished as 
follows: Express the last term as a sum of two terms in the order in which 
they appear; then from each take out a factor corresponding to n” = n 1 and n" = n 
respectively. After the indicated factor is removed, the first factor becomes 


i 

3 


E E 

n'^n n"^ £ n,n' 


H^,, H ( t> + HO), H<*) , 

n ' n 1 < rv»rt nrv» n'm ' n 'n 


n'n 


n”n 


nn M 


n M n' 


n'n 


(E° - E° J (E° 
n n T n 


E° ) 
n 



E 

n'#n 


(H (l) , H^) H^, 

nn 1 n'n n'n' 

(E° - E 0 ,) 2 
n n' 


and the second factor becomes 


a 1 E E 

n'^n n"^n, n' 


nn T n'n' 1 


H ( ^ +H (1) „ 

n"n nn" 


H 


(i) 


n"n T n T n 


h <:>) 


(E° - E 0 ,) (E° - E° ) 
n n ?/ n' n"' 


2 

3 


E 

n^n 


nn f n r n' nn 


o o 9 
(E - E ,) 2 
n n' 


Note that the last two terms of the above expressions combine with the third 
term of equation (30). Thus, equation (30) is modified as follows* : 


e (3) _ H ( 3 ) + 


n 


nn 



(H (2) , + H (1) , H ( f ) H (1) , H ( [ ) (H ( J\ - H (l) ), 

nn r n'n nn T n'n . nn T n T n n'n' nn 

■ *f* 


(E° - E°,) 
" n’ ' 


(E 


n 


E °,) 2 

n' 


* This result agrees with equation (11-33) in Powell and Crasemann's book on 
Quantum Mechanics. 
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E 

n'^n 


E 


(H (1) , 

nn' 


H 


(i) tt( 1 ) . h (D H (i) 

n'n" H n"n + H nn" H 


n"n' 


H^) 

n'n 


n'^n, n' 


(E° - E 0 ,) 
' n n ' 7 


X 



( 31 ) 


A more compact form is obtained by expressing the double summations in a dif- 
ferent manner. Take the above term and write it as two terms after inter- 
changing n" and n T in one of them. This gives 


1 

3 


E E 

n'^n n"#n,n' 


H (1) , H^,, 
nn 1 n'n" n"n 


(E° - E°,) 
' n n ' 7 


0 0 0 „0 
E - E „ E - E 
n n" n' n" 


£ E E 


H (l) , H ( J\, H ( ^ 
nn' n'n" n"n 


a u lj o o 

n'#n n"^n,n' (E - E^,,) 


o o 

E - E , 
n n' 


o o 
E . - E 


n' 


n' 


Combining these two terms again one gets 


1 y y „(•)„(.) H (.) 

3 Jin n"*n,n' n ' n " n " n \<E° - E° ) (E° - E° ) (E* - E°„) 


tr(0 „(*) H (l) 

ri , Jtl , n ,, 
nn' n'n" n"n 


- E E 

n ’* n n " #n ' n <-E^,) <E°-E»„) 


Now, one puts back the factor corresponding to n" — n' in the above expression 
and obtains 


H (i) w (i) H ( 1 ) 

nn' n'n" n"n 


(H (1) , H^) 

' nn' n'n 7 n'n' 


Z y mr_ in i' ' n ii y 

Zj nnoo Li OOo 

n'^n n"*n (E n - E^) (E n - E n , f ) n'*n (E n - E ^) 2 


I II II ■■■ II I III I II 


i ip ■■■phi ■Him hi i i mi u 


i ai 


When this result is substituted in equation (31) , the following is obtained: 


e (3) _ H ( 3 ) + 

n nn 


E nn 1 n'n nn' n'n _ ^( 1 ) y 
n n nn Li 


H (o H ( ;> 

nn' n'n 


n'^n 


(E° - E° ) 
v n n' 7 


o o o 
n'^n (E - E ,) 2 
n n' 7 


nn' n'n" n"n 


7 7 

n'^n n"#n (E^ - E^,)(E^ - E® It ) 


(32) 


The above result is easier to use in actual calculations. 


V. FOURTH ORDER FORMULAS 

( 4 ) 

In this section, the expression for E^ is derived. From the general 
form developed in Section II, one obtains the fourth order correction to the 

energy. As was the case for 7 , one needs to substitute known quantities foi 

[S^7 H^] , [S (2) , H^] , and [S^, H^]. Some of these quantities were 
defined in the preceding section. From equation ( 27) one solves for the re- 
maining expression: 


[S 


(3) ; H (0) ] = (E ~ ~- H ^ ^ ([S (2) , H (l) ] + [S (1) , H (2) ] + [S (2) , E (l) ] 


+ [S (1) , E (2) ]) - ~ ([ S (l) , [S (1) , H (l) ]] 

- [S (l) , [ S (1) , E (l) ]]). 


(33) 


In equation (21) note that the term above is also found within brackets. Hence, 

in addition to the above, one needs to calculate [S^ \ [S^, H^]] . The ex- 
pression for this term is 

- [S (I >. S^-] [S (!) , H (,) n +[S (1) , [ s <1 > ,h( 2 >] ] 

1 1 Ci 

+ [S (l) , [ S (2) , E (1) ]] +[S (1) , [S (l) , E (2) ]]) 
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- (tS (l) ,[S (l) ,[S (l) , H (l) ] ] ] - [S (1) ,[S (1) ,[S (1) , E (1) ]]])} . 

With the above expressions substituted into equation (21), E ^ becomes 
E (4) = H (4) +i{[S (4) ,H (0) ] + [S (3) H (l) ] + [S (2) , H (2) ] + [S (l) , H (3) ]} 

+ |[S (1) ,E (3) ]-|[S (1) ,H (3) ] +~ [ S (l) [ S (2) , H (l) ] ] +“[ S (l) , [ S ( t h <2) ] ] 

+ ^[S (1) ,[S (2) ,E (1) ]] +|[S (1) ,[S (1) ,E (2) ]] +^[S (1 >[S (1) , [ S (l) , h (1) ]]] 

-■^tS (1) ,[S (1) ,[S (l) , E (1) ]]] -^[S 2 , +|[S (2) ^~] 

+ ^[S (2) , [S (l) , H (l) ]] +^[S (2) , [S (l) , E (l) ] ] -|[S (3) ,^]+f[S (3) ,^] 

-^[S (1) ,[S (2) , H (l) ] ] - |[S (2) ,[S (1) ,H (1) ]) -|[S (1) , [S (l) , H (2) ]] 

" 6 t S (1> [ S^ l) ,^T~] ] + |[S (1) ,[S (1) , fi |p] +^[S (1) ,[S (1) ,[S (1) ,H (1) ]]] 


+ T 2 rs (1) , [s (1) , E<‘>,„ [S < 2 >, 

-^[S ( 1 ) ,[S ( 1 ) ,[S ( 1 ) ,H (1) ]]] -^[S ( 2 ) ,[S (1) ,^-]] + j[S ( 2 ) ,[S (l) , ii ^-]] 

+ ^[S (l) ,[S (l) ,[S (l) ,^-]]] -^[S ( 1 ) ,[S ( 1 ) ,[S (1) , ^]]] . 

Upon recombination, this becomes 


r* 1 ) 


E^ 4 ) = H (4) +|{2[S (4) , H (0) ] + [S (3) ,H (1) ] + [S (2) , H^ 2) ] + [S^, H (3) ] 


+ [S (1) , E (3) ] + [S* 2) , e ( 2) ] +[S (3) , E (l) ]} 
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( 34 ) 


-£{rs ( 1 >.[s (1 >, 


H (l) ]] + [S (l) ; 


[S 


(i) 


H (2) ]] 


- [S 


(l >, [S (2 > 


E (1) ]] 


-[S 


(1) , [S (l) 


e (2) ]] 


f [S (2 >, [S (l) ; 




-[s ( 2 ) ,[s (1 ) ; 


E (1) ] ] } 


This result is quite an improvement over the previous expansion. 

(4) 

Once E^ is expressed in its simplest form, the diagonal terms are cal- 
culated. Expanding the brackets, one gets the following terms for the diagonal 
(excluding the zero elements): 


£< 4 ) =H (4 > + f 
n nn 2 


(3) , H (l) ] + [S (2) ; 


H (2) ] + [S (1) : 


h (3) ]} 


nn 


+ J-{(S (1) h ( 2) S ( 1 ) +s (l) h ( 2) S ( 1 ) ) _ ( H < 2 ) S (l ) g^ 1 ) + g(l) g^ 1 ) H ( 2 )) 

X iLi 

+ (S^ E^ - E^ S^) + (E^ E^ S^) 

+ (S (2) H (1) S (l) +S (1) H (1) S (2) ) + (S (2) H (1) S (l) +S (1) H (l) S (2) ) 

- (H (1) S (l) S (2) +S (2) S (1) H (1) ) - (S (1) S (2) H (1) +H (1) S (2) S (l) ) 

+ (S^ S (l > E^ + E^ S (1) S^) + (E^ E^ ) 

- (S (1) E (l) S (2) + S (2 > E (1) S (1) ) - (S (2) E (1) S (1) + S (1) E (l) S (2) )} 

* nn 

(35) 

(3) (4\ 

Except for S , , all the terms necessary to calculate E 7 have been 

n'n , n 

(3) 

calculated previously. Before one can continue, must be obtained. It is 

obtained from the nondiagonal terms of equation (33). These terms, obtained by 
expanding equation (33), are: 
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[S (3) , H(0) ] n , n = i -|{[S (2) , H (l) ] + [ S (l) , H (2) ] +fS (2) , E (l) ] 


+ [S (l) , E (2) ]} + ~ { (S^ H (l) S (l) - S (l) S (l) H (l) ) 

1 J J n T n 12 1 ' 1 


+ (S^ S^) + (s^ 1 ^ E^ 

- S (l) E (l) S (l) ) 

+ (E (1) S (l) S (1) - S (l) E (l) S (1) )} n , n . (36) 


!<?) 


The above equation is arranged so as to simplify the calculation of S^t^. This 
fact is realized only after trying various combinations. Since 


[S (3) , H (0) ] , = Yj ( s( ?\, h( m - h( ! ) „S ( ?, ) ) 
1 J n'n t-i n'n" n"n n'n" n"n 

n" 


= ( E° - E° ) S ( ?> , 
n n T/ n'n 

one can solve for S^ 3 ^ from equation (36) by simply dividing by ( E^ - E^ t ) . 
This step will be carried out after each term in equation (36) is evaluated. 

Before proceeding with the calculation of , one should modify S^ 
v n'n J n'n 

to simplify matters later on. To do this, one separates the third term of 

equation (28) into two terms and redefines the resulting terms as follows: 


E 


H ( ?>. . H (1 > 


E 


ii 0) „<»> 

n'n" n"n 


H*}*, H ( ^ 
n'n* n'n 

n'^n (E° - E° ) (E° - E° n ) ln"*n,n’ (E° - E°.) (E° - E°„) (E° - E°,) 2 

n n 1 ' ' n n"' n n" n n"' n n 1 ' 


E 


„(i) „d) 

n f . xi 
n'n" n"n 


E 


„<1> H (l > 

Jtl # .. XI tf 

n’n" n”n 


H ( ;> h"' 

n’n nn 


n"#n' (E - E ,) (E - E „) ln"^n,n' (E -E,)(E -E„) (E - E ,) 

n n f n* r» M/ v ' n n T ' n n n/ ' n n * ' 


n T n 


n 


n 


n n 1 


One can see that the factors removed from the sums combine with the second 
term in equation (28). Incorporating these changes, becomes 
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i H 


( 2 ) 

n'n 


3 (!) = 

n ' n (E'-E 0 ,) 
' n n' 7 


i H 


d) /„<*> 


( 1 ). 


n’n 


(H , - H ') 

n'n' "" 7 


nn 


(E° - E°,) 2 
n n ' 7 


5 Z 


H 


(1) H (i) 


n'n" 


n"n 


O 1_J 0 0 

n"*n,n' (E - E ) 
n n' 7 




n 


n' 


E° - E c 
n' 


n' 


( 37 ) 


( 3 ) ( 2 ) 

In Section IV, E was obtained without making this modification in S , . The 
n & n'n 

increased complexity of the present case makes it imperative that this change 

be made here. 


With the modified form of , the calculation of S , proceeds smoothly. 

n'n n'n J 

In the calculation of the terms of equation (36) , one will find it more convenient 

to modify each term as soon as it is calculated. If this is not done, the situation 

becomes complicated. Taking the second term of equation (36) one obtains the 

following: 


1 (2) (1) _1 y (l) s ( 2 ) _ s (2) (1) 

2 [ ’ 1 n'n 2 Lj (H n 'n" n"n n'n" n"n ) 

n" 


„(1) TT ( 2 ) 

i ; V n'n" n"n 

~ 2 \ Z + 


z 


„(») „(1> 

n'n" n"n 


2 ^n"^n (E° - E° ) n"#n' (E°, - E°„) 

v n n" 7 n' n" 7 


Hd) nd> (hW -Hd)) nd) nd) md> -nd) ) 

^ n'n" n"n v n"n" nn 1 + ^ n'n" n"n K n"n" n'n 1 ' 

^ ( E° - E° ) 2 n"#n' ‘ * 2 

' n n" 


n"^ 


{E n’- E n»>‘ 


ill l 


(i) (i) (l) 

n'n" n"n'" n"'n 


i"*n n'"*n,n" (E° - E°„) 

n n" 7 


oo 0 0 

* E n“ E n'"* (E n" “ E n' n) 


z z 


( 1 ) ( 1 ) ( 1 ) 

n'n’" n"'n" n"n 


n"*n' n'"#n\n" (E° -E°„) 

n' n" 7 


. 0 0 
\E . - E 


0 0 
E - E 


n' n"' n" n' Mi 
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By taking out appropriate terms from the above sums, one can write these 
terms in a more symmetric form. The steps necessary for this modification 
are as follows: Take out the factors corresponding to n" = n' in the first term 
and the one for n" = n in the second term, this gives 


I 2 


nf 1 ) tt( 2 ) jj( 2 ) ) 

n'n 11 n"n n'n" n"n i n'n 1 n'n-' nn ' 

— + I • ,a) 


n"#n,n' \ (E n - E n „) (E nl - E n „)/ (E n - E nt ) 


Repeating this procedure with the next two terms one obtains 

jC 1 ) tt(*) ) 

n'n" n !! n ^ n n n M nn ' n'n" n M n ' n"n" n'n'* 


i 2 


n'^n.n' 


(E° - E 0 ,,) 2 
n n" 


< E n' ' E n "> 2 


. h"* (H ( I\ - H (1) ) 2 

1 n' n n'n* nn 

2 0 o 2 

(E - E ,) 

' n n* 7 


(b) 


If this step is repeated with the last two terms, the result is 


1 J 

f H d) H (l) H (l) 

1 y ^ n'n" n"n'" n"'n 

l i 

. i \ 

4 } 

, h . L o o 

i 0 n 0 _o 

\ 

[n"#n,n' n"'*n,n" (E “ E n ..) 

' E n - n"' 

n" n"' 1 


2 2 


(i) (i) (i) 

n'n'" n"'n" n"n 


n'^n, n' n'"#n',n" (E^, - E^„) 


E° - E , 

n' n' 


E°„ - E °i 
n" n' 


I 2 


H< n (h^.-h* 1 *) 
n'n" n"n n'n' nn 


n"#n,n' 


o o 

(E - E ,) 
n n' 7 


o o 

E - E 
n n" 


o o 

E - E 
n* n' 


(c) 


After the individual terms are obtained by removing terms corresponding to 
n" = n' and n" = n from the main double sums, the last term in (c) is obtained 
by letting n'" = n". Thus the second term in equation (36) yields terms (a), 
(b) , and (c) . These terms will be incorporated with the other terms later on. 
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The remaining terms of equation (36) are given in their final form only. It is 
assumed that the reader can retrace the source of the extra terms generated. 
These terms are: 


J[S (1) , h (2) ] 


/h ( ; ) „h ( ;, ) h ( ;\,h ( ?m h, (h i :\-h' 2 ' 

_ _i y I n'n" n"n f n'n" n"n I { 1 n'n n'n 1 nn 

n '” = 2 n4„,„. (E* - E“„, <E°, - E°„)/ 2 (E° - E \) 

1 v n n" 7 ' n' n" 71 ' n n' 7 


( 2 ) „( 2 ) 


(d) 


1 , . . (E^ - E (1) ) . (E ( ^ - E (1) ) 2 

1 g ( 2 ) E (i) i n'n n' n ’ i n’n n 1 n 7 

2 ’ n'n 2 (E° - E° ) 2 

v n n' 7 
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n n' 


i Z 
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n'n" n"n n'n' nn 


n'^n, n 1 


(E° - E 0 ,) 
n n' 


1 1 
+ 


\E° - E° E u - E w 
n n" n' n" 


» (e) 


-rs (1) e (2) i - 

2 [S ’ E n'n ~~ 2 


(E ( p - E (2) ) 
i n' n n'n 


(E° - E° ( ) 
n n' 7 


(f) 


~ (S (l) H (l) S (l) - S (l) S (l) H (l) ) 

n'n 


-I z 


Til') „(1) II (!) 

rl il 

n'n nn TI n M n 


6 n'^n (E°-E° ) (E°-E° H ) 12 

' n n' 7 v n n" 7 


i Z 


„(D H (l) m (l) H (i). 

n , n (Jti - ri ) 

n T n n n M n v n M n" nn 


n'^n'.n (E° I -E°„) (E°-E° M ) 

n' n" v n n" 


k Z Z 


(i) (») (i) 

n'n" n"n'" n'"n 


12 ^ t- 1 0 0 

n'^^n' n'"^n,n" (E , - E ,,) 

n' n" 


1 1 

+ ■ 


F° F° f° F° 

E n E n"' E n" E n"" 


(g) 


y- (S (1) H (l) S (1) - H (l) S (l) S (1) ) , = 
12 ' 7 n'n 
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H 


(1) H (J) „(1) 


in" 


n'n 


n M n' 


H' 


n'n 


+ i v -- — — -t- 

6 n"*n’ (E°-E°,) (E° -E°„) 12 

n n'' ' n f n u/ 


Z 


„(*) H (i) m (i) H (l) v 
H . .. H .. (H n „ n „ - H n , n ,) 


n'n' 


n"n 


n"*n,n' (E° - E° ) (E° - E° ) 

' n n" 7 n' n" 7 


t 2 l Z 

n'^n, n* ^"T^n'in" 


(l) (l) (1) 

n’n'" n'"n" n"n 

(E° - E° ) 
n n" 


K - <•" 


n' 


- E 


n’ 


(h) 


_L , s (*) c (l) 

12 1 


S (1) E (1 > -S (l > E (1) S (1) ) Y 

; n'n 12 . fJ 


H 


(i) „(!) ^ 

n'n" n"n ' nn n"n" ; 


n"*n,n' (E° - E°„) (E°, - E°„) 

' n n" 7 ' n’ n" 7 

(i) 


~ (E (1) S (l) S (1) - S (1) E (1) S (1) ) 

12 'n f n 


i 

12 


i 


H^,, (H^, - „) 

n'n" n"n n'n' n"n" 7 


n^n,n' (E° - E°„) (E° , - E° „) 
' n n" 7 v n 1 n" 7 


( 3 ) 


Combining (a) , (b) , and (c) with all of the terms just obtained above and 

dividing equation (36) by (E° - E°,)/i, one obtains the following for s(?): 

n n’ n f n 

(cont'd on next page) 
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3 (3) H (3) H (3) (H^, - H (1) ) 

n'n n'n n'n v n'n' nn 

i ~ o 0 0 _,0 2 

< E n - E n’> <E n ' V 


1 H n’n H n'n’ im ’ + 2 
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/H^ 
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' n n’ ' ' 


n'n" n"n n'n 1 nn 
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(E° - E° ) 

' n r» T ' 


H (l) rr(i) H (l) 

y n'n" n"n'" n"’n 
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H (i) H (i) H (D 
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+ L< o o o o 

n"'*n',n" 4(E -E ,) (E ,-E „) 
' n n' ' n' n" 


(1) (1) (1) 

y n'n" n"n'" n'"n 

Li n o o o 

n"'*n,n" 12(E°-E U ,)(E ,-E ,,) 
' n n" n' n" 


ttO) «(*) fj(i) 

+ xp n'n'" n'"n" n"n 

n»'#n»,n" 12(E°-E 0 ,) (E°-E°„) 
’ ' n n' 7 v n n"' 
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A few comments are necessary in regard to the above result. Terms (i) and 
(j) were cancelled by the middle terms of (g) and (h) . By expressing the first 

/2) 

terms of (g) and (h) in terms of 7 , one can combine them with the second 

term in (d) and the single term in (f). The last modification is accomplished 
by taking the first terms in (g) and (h) and expressing them as follows: 


By taking 




n"#n' (E°, 
n' 


h'hi* 

E° ) 
n"' 


/ H 


( 1 ) 

n'n 


< - E l 


_i 

6 


H 


(l) 

n'n 


tt( 1 ) „(*) 

11 n 
nn" n"n 


0 0 t-J 0 0 

(E -E ,) n"*n (E u - E u „) 
v n n' 7 v n n" 


then using 


E ( 2 ) - jj^ 2 ^ + 
n nn 


l 


„(i) „(i) 

n , , n 
nn" n"n 


n'^n (E° - E° ) 
' n n" 7 


one obtains 


_i 

6 


H 


^ (E^ - E^) 
' n ! " 


n*n 


n 


o o 

(E - E ,) 
v n n' 7 


_i 

6 


H 


(1) ,„(*) 


( 2 ), 


n'n 


(H\ 7 , - H' 7 ) 
n'n' nr * 


nn 


(E° - E 0 ,) 
v n n' 7 


These terms combine with (f) and the second term in (d) , respectively. 

With the expression for established, equation (35) is now evaluated. 

n'n 

The same procedure is followed; that is, each term is evaluated separately and 

then combined to form E ^ . Going back to equation (35) , one notes that the 

first term to be evaluated is f , H^l . Using the expression just obtained 

1 nn 

( 3 ) 

for S , , one obtains the following: 

n'n 
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(H (1) , + H (1) , H ( f ) 

Im , n’n" n n n nn T n f n u n n n 

4 (E° - E° f ) 

' n n' 7 


E° - E°. 


E° - E° 
n' n' 


mf 1 ) w (2) xt i 1 ) . tt( 2) H (l) H ( 1 ), 

^ nn n n T, n T n ! n nn n n n n T n'n 
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+ 


z 


n"' ^n, n" 


(l) (i) (i) (i) 

rm" 1 n"'n" n"n , n'n 

24(E° - E 0 ,) (E°, - E° ) 
' n n' n' n" 



The range on n'" is not the same on all of the above terms. For this reason 
the terms containing n'" are not fully combined at this point. Without making 
additional simplifications, it is helpful to rearrange the above term as follows: 
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(4) 

The above terra is the longest and most important in E ' . The remaining 

terms are obtained as a group below. Omitting intermediate steps, these terms 
are: 
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n'n'" 

n'" 

0 

n 

E - 
n 

E ») 
n' 7 


0 x. 0 

E - E , 


o 

E - E 


n' 


U J 


(0) 


■ jr (S (1) S (2) H (1) + H (1) S (2) S (l) ) = 

12 ' "" 


nn 


2 

n'^n 


(H (2 (h ( ? +h<1 ! h( ?)H (1) H (1) , H ( 7 (H ( J\ - H (1) ) h (1) 

1 ' nn* n’n nn* n'n 7 nn 2 nn' n' n n'n' nn nn 


12 


(E° - E°,) 2 
n n' 7 


12 


( E° - E 0 ,) 3 
n n' 7 


- 2 
n"*n,n' 


H (3) , H ( 7 +H (1) , H ( ^„ U^} ) 
1 nn” n"n' n'n nn' n'n" n"n 7 


12 


(E° - E°,) ( E° - E° ) 
' n n' 7 n' n” 7 


(H (1) | H ( | ) „H ( J ' 

1 nn' n'n" n"n nn" n"n' n'n 7 n"n" n'n' 7 


12 


0 0 0 0 9 
(E n' E n’> < E - 


n' 


/„(*) jt (1 ) H ( J ) . H (i) H (i) H 0)v „(!) 
1 nn' n'n" n"n nn" n"n' n'n nn 


24 


(E° - E°,) 2 
n n' 7 


x , 0 xx 0 

E - E 
n' n' 


E ° - E °„ , 

n n" / 


( H nn* H n'n'" n'"n" H n''n_~n n" _n"n’" n"'n' n'n 


- 2 

!4 n"'*n',n" (E^ - E^,) <E», - E;„) 


0 "0 x /P o ~0 


0 0 0 0 
\ E n' - E n"' E n" " E n'" 


, (0 
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~T (S (2) S (1) E (l) + E (l) S (1) S (2) ) 

12 'nn 


Z 

n'^n 


u 

(h (2) , 

nn' 

H^J) + h^ 1 ) H < 
n T n nn f r 

r 2) ) H (1) 
m nn 1 

L 12 


0 0 2 
(E n - V 

+ 6 



vi nn" 

H ( ,V , H ( P 

n"n' n'n 


24 

L 

n'^n, n' 

(E° - 
n 


(E° - E°,) 3 
' n n" 




0 0 0 

E - E E - 
n n" n' n 


J] 


(k) 


jr (E (1) S (2) S (1) + S (1) S (2 ) E (l) ) = 

12 7 nn 


Z 

n'*n 


1 nn f n'n nn' n'n 7 nn 


12 


(E - E ,)‘ 
n n' 7 


H (1) , (H*J\ - H (1) ) H (1) 

1 nn n'n n ( n' nn 7 nn 

6 ( E° - E 0 ,) 3 

n n' 7 


— 7 

24 „p 


(H (l) , H ( j\ t H ( y +H (1) n H ( y H ( J } ) H (1) 
nn 1 n'n" n"n nn" n"n' n'n ; nn 


n'^n, n' 


(E°-E° t ) 2 
n n' 


1 




(\> 


(The last two terms yield the same results.) The next two terms are identical; 

E^ 1 ^ + S^ E^ S^) 

12 7 nn 


Z - 

n'^n L 


(2) „d) 


H\' + H (l) . H (2) ) H ( !\ 
2' nn’ n'n nn' n'n 7 n'n' 


12 


(E° - E° ) 2 
n n' 


H (1) , /H^, - H (1) )H ( *\ 

2 nn' n'n n'n' nn 7 n'n' 

6 <e°-e\) s 

n n' 7 
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m ( 1 ) H (1 > H (l > + H (l > H (l) H (l) * 
2 v' 1 nn" n"n* n’n nn’ n'n" n"n' 


n"^n, n' 


( E° - E°,) 2 
n n' 


1 1 

+ ~ 

E° - E° E° - E° 
n n M n 1 n M 


To simplify the algebra, the above terms just obtained (a, /3,y, . . . ) are com- 
bined separately before they are added to H^ 4 ) and [S^ 3 ) , H^ 1 ^] . Taking terms 
a, /3, y, . . . one obtains the following expression: 


s r ( 


+ H (3 >, H ( :>) 

nn' n'n tin 1 n'n' nn' n'n 

2 (E° - E° ) (E° - E°,) 

' n n" n n !/ 


H (1) , (H ( J\-H (2) ) hW,H<? (E ( J } -E (2) ) 

nn T n T n n T n T nn ' nn ? n T n n T n 


6 (E° - E 0 ,) 2 
x n n* 


6 ( E° - E° ) 2 
' n n 1 


(H (1) , (H ( |\ - H (l) ) 

mr n T n nn T jrn n ? n T nn 

(E° - E° ) 2 
n n* ; 


Z 

n n ?%i, n T 


(H (l) „ , H (2) + H (2) , H ( ^„ H ( ^ ) 

nn M n"n T n In nnl n T n n n n n 

3 (E° - E° ) 
n n* ' 


1 1 

+ — 

E° - E 0 „ E° - E° 
n n M n 1 n M 


(H (1) , H (2) „ H<J> + H (l) „ H (2) , H^) 
nn T n T n n n Tt n nn" n M n T n T n 7 

12 (E° - E° t ) 
n n T 


E° - E° E° - E° 
n n n n 1 n f 


(H (2) tJ H ( ^ , + H (1) , H^,, H (2) ) 

n n IJ n’ 1 !! 1 n’n nn 1 n'n” n"n 

o o 

12 (E - E ,) 


E° - E° E° - E° 
n n" n' n" 


(H (I) , H^„ H ( ^ + H (l) „ , H^) ( H ^ni " H !f ) 

nn ! n T n ,T n ,T n nn n n tf n T n T n ; n n nn 

6 (E° - E 0 ,) 2 (E° - E° ) 
n n' n n" 


I 



(H (1) , H*)\, + H (1) n , H^) (H ( j\ - H (l) ) 

rm T n f n M n n n nn Tt n"n T n ! ir n T n T nn 


0 0 9 0 0 

12 < E „ - V < E m ‘ E »"> 


ir n T 


(H (1) , H^,, H ( ^ (H^V - H ( }\) 

nn* n'n" n"n nn" n"n , n'n n"n" n'n' 

oo o 0 2 

12 (E n- E n’> < E n'" E n"> 

2 +H (1) , H*J\ f H ( ^ ) 

v nn IM n"'n" n"n' n'n nn * n'n" n"n'" n'"ir 

24 (E° - E°,) (E° - E° ) 
n n' n n" 


* 2 
n'"^n,n" 


/ — - — + — 

1e° - E° E° - E° 

\ n n m n" n' 


♦ 2 

n' M #n, n' 


m (I > H (l) H (l) H (1) H (1) H (1) H (1 > H (l > t 

nn 11 n"n t n'n"' n ,M n nn"' n" , n , "n'n" n"ir 

24 - E J"> ( E n'- E n'.> 

X ( i 1 

\E° - I 


-t- ■ 


E , 
n n' 


E°, - E° 
n' n' 


* 2 

n'" : » t n',n" 


, (l) (l) (l) (l) „(l) „(l) „(l) H (l) 

nn f .n , n ,n n fM n n n n n nn n n M n™ n MT n T n T n 

24 (E° - E 0 ,) (E° - E 0 ,,) 
v n n' n' n" 


1 1 

+ 


e° - E® 

E n' “ E n'" n" n m 


Incorporating the above terms with [ S ^ , II ^ ^ ] and one obtains, 
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= 

n 


h (1) + E 

nn u 


h (2) , h (2) h (1) , h (3) + h< 3 > h<!> 

nn 1 n ! n nn’ n f n nn' n'n 


ft’*n \ (E» - E» ,) 


0 0 

(E - E ,) 
' n n' 7 


nn’ n'n ' n’n' nn ’ v nn' n'n nn* n'n 7 . n'n* nn 


(E° - E°,) 3 
' n n' 7 


(E° - E° ) 2 
' n n' 7 


H “> h<:> (E <:>-e< 2 >) „<•> h<:> ( h ( ?>,-h< 2 »,- 

nn ' n'n ' n' n 7 nn' n'n n'n' nn 


2 (E° - E° ) 2 
v n n 17 


0 \ 2 


2 ( E° - E° ) 
n n' 


*2 S , 

n'#n n’^n, n’ 


H^„ H (2) + H (2) „ , H^) 

' nn T n ! n" n n n nn ,T n M n T n'n 


3 (E° - E®,) 
n n' 


o 

E - E 
V n n" 


E° - E° 


n" / 


( 1 ) ( 1 ) ( 2 ) ,„( 2 ) ( 1 ) ( 1 ) 

H nn” n" n' H n'n H nn' n'n" H n"n ) 

3 (E° - E° ) 
n n' 


I 1 

E° - E° 

^ n n" 


E , - E l 
n' 


n' 


H ( ?„ H ( J +H (1) II H ( J | H ( [ ) ) / 

nn' n'n" n"n nn" n"n' n'n 7 


o o 

3 (E - E ,) 
n n' 7 


E° - E°„ E°, - E° 

\ n n" n* n” 


nn' n'n" n"n nn" n"n' n'n 7 n'n' nn 7 


nn 


3 (E° - E^ f ) 2 (E® - E®„) 

nn T n t n tT n M n nn M n !T n' n'n 7 x n'n' nn 7 
3 (E® - E^,) 2 (E®, - E®„) 

(H (1) , H^,, H ( ^ +H (1) „H ( J, ) H^) (H ( ^ „ - H ( J\) 

nn' n'n" n"n nn" n"n' n'n 7 n"n" n'n' 

3 (E® - E 0 ,) (E° - E 0 ,,) 2 

n n' 7 n' n" 
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e 

n'"#n, 


E 

n"'^n', 


E 


E 

n'"^n', 


E 

n'"#n, 


) 

5( nn' H n'n" n^n"' n'”n nn’" n' "n" n"n' n'n' 


24(E»- E ',) (<-e;„) 




E 

n n’ 


o i o 

E - E 
n" n' 


^ /tt (l) Tj(i) tr(^) tj(*) + Jj(*) JJ^) ) 

4 E nn" **n"n'" E n"'n' E n'n nn* n'n'" n" l n" n"n_ 


n’ 


24 (<-<•) < E n." E n») 


1 1 
+ - 




,„(*> H <‘> H (1) H (1) +H (1) H (1) H (1) H (1) > 

( nn' n'n" n"n'" n'"n_ nn'" n’^n" n"n^ n’n ; 


24 ( E° - E°,) ( E° - E° ) 
n n' n' n" 


On oo 
E n " E n’" E n" ” E n"" 


CH^ 1 ^ H^+H^H^ H (l) H (l) ) 

nn" H n"n ,M n'"n' n'n + nn' n'n'" _ n'"n" n"n' 


n' 


24 (E° - E° ) (E° - E®) 
' n n T n n n/ 


0 0 o 0 

> E n'-V E "- E - 


H (1) H (,) H ( ‘> +H (1) H (l> H (1) H (1) 

1 nn" n"n' n'n'" n"'n im'" n"'n’ n’n" n"n 


I 


24 (E° - E° ) (E° - E° ) 
v n n ' 7 n' n"' 


L E° - E E° - E° 
n n'" n' n 


HI 
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If the above result is to be used for numerical calculations, it is neces- 

(4) 

sary to simplify E as much as possible. With the exception of the terms 

(4) 

containing triple sums, all the terms in E^ can be recombined into a more 

compact form. The range on n"' may be standarized by removing terms 
corresponding to n' M = n, n', or n" where necessary in order to obtain quantities 
of the form 

2 2 Z «— • 

n'#n n"^n,n' n" ,: #n, n', n" 

(4) 

Referring to E^ 7 , one needs to take out the factor corresponding to n"' = n' 
in the first triple sum which gives 


_5 

24 


Z Z 

n'#n n'^n.n 


(H ( ^, 

nh' n'n' 


H ( :> , h^) 

n"n' n'n 7 


(E u - E ,) 
n n' 7 


(E 


n 




n’ 


.) 


E u - E" 

V n 


n' 



(A) 


From the second term, one needs to take out the factor corresponding to 
n m = n, which gives 

(H (1) „ H ( * } H (1) , + H (1) , H (l) „ ) 

_4 y nn" n"n nn' n'n nn 1 n'n nn" n"n 7 

24 4 1 P 0 „0 . ,_,o n 


n'^n n'^n, n' 


(E - E ,) (E , - E ) 
n n' 7 n' n" 7 


0 0 

E , - E 
n' n 


E - E 
n" tf 


Letting n" = n' in the second factor, the above reduces to 


tf Z Z 


nn" n"n nn' n'n 


24 n'^n n'^n'.n' (E° -E°„)(E° - E° t ) \ E u - E° , E° - E°„ 

' n n" 7 n n' 7 ' n n* n n" 


_i 

r»0 _ xt'O 


(B) 


From the third term, for n 1 " = n\ one gets 


Til Z 


,„(*) H (l ) H (l > H (l >\ 
nn' n'n" n"n' n'n^ 


12 lj o o o o 

n ! ^n n”^n,n ! (E -E f )(E° -E „) 
x n n T/ ' n f n"' 


E u - E l 
n r 


E° - E° 


n' 


n" / 


(C) 
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From the fourth term, for n' M = n, one obtains 


1 y V 1111 n"n 1111 n'n 

12 4 / Lj 


n’*n n"*n,n» (E -E „)(E - E ,) \E° -E°, E° - E° 

' n n" ' n n ,; \ n n* n n"/ 


(D) 


From the last term, for n m = n", one obtains 


nn' n'n' 1 n"n' n'n 7 

* 2 n'^n n"#n,n' (E° - E° )(E° - E° ) 

' n n" 7 n' n" 7 


-s Z .. Z 


E° - E°, 
n n 1 


E°. - E 


n* 


n" / 


(E) 


Combining the above terms, one obtains the following: 
From (C) and (E), one gets 


(i) (i) (1) (1) 

im’ n’n" n"n' n' n 

12 n'*n n"*n^n' (E° -E° f )(E° - E°„) 
n n' n n" 7 


n S .1 


E° - E° , E° - E° 

. n n' n' n 


This, in turn, adds to (A) giving 


1 y y nn 1 n'n" n"n' n'n 

2 n'^n n"*n,n’ (E° - E°,) (E° - E°„) Ie° - E°, 

n n' 7 ' n n" 7 \ n n’ 


E° . - E° 


n‘ 


n’ 


(ACE) 


Terms (B) and (D) yield 


-iS S 


(H (1) „ H (i) , H ( J } ) 

nn" n"n nn' n’n 7 


n’^n n"*n,n> (E^ - E® ,) (E* - E* lf ) j E n - E n , E^ - E« If | 

.( 4 ) 


(BD) 


Terms (ACE) and (BD) must be added to E 7 when it is expressed in terms of 
the newly defined triple sums. 


Since the triple sums are now alike in their summation indices, one 
recombines these terms into a more symmetric form. The most logical arrange- 
ment would be to express each one of the terms as matrix products of the form 
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H (l >. H (1 > 


H 


(i) H d) 


nn' n'n" n"n’" n'"n 


In order to accomplish this, consider the terms in the order in which they ap- 
pear in E^ . Multiplying the first term out and separating the product into 
four factors, one obtains 


/H (1) , H*}*,, H ( J f 

_5_ y y y / nn n*n" n"n"* n' M n 

24 Li Li Li \ oooooo 

n’*n n’^n.n' n"'*n,n\n" l(E -E ,)(E -E „)(E -E „,) 

» ’ ’ *s' n n' 7 n n" 7 n n"’ 7 


(1) (1) (l) (i) 

nn’ n’n" n"n m n"’n 


(E° - E° ) (E° - E° ) (E° - E 0 ,,,) 

' n n’ n n" 7 ' n" n M,/ 


n' 


n n' 


n' 


H (i) H (i) H (i) „(l) 

nn'" n n, n" n"n' n'n 


0 0 0 0 0 0 

(E - E ,) (E - E ) (E - E ) 
n n' 7 n n" 7 n n"' 


(i) (1) (i) (i) 

iy ... n ... .. n .. . n . 


(E° - E° ) (E° - E° ) (E° - E° ) 

n n* n n" 7 n" n"' 7 


The first two factors appear in the desired form, the last two do not. In order 
to modify them, one interchanges n'" and n’, and obtains 


V 1 ) H (.i) H (l) H (i) 

nn’ n'n" n"n'" n"’n 


24 Z; \ oo oooo 

n'^n n M ^n,n' n" , #n,n , ,n" \ (E -E ) (E -E )(E -E ,) 
> ’ ’v' n n"' 7 n n" 7 n n' 7 


H (1) . H ( |> H (l) 


nn T n T n u n M ir 


H 


(i) 


fn" n n n nT n"» 


n MI n 
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Applying the same analysis to the rest of the terms, one obtains the following 
set of factors* (including the ones above) : 


z z z 

n'^n n"^n’,n n' M ^n",n',n 


10 


|_24 (E° - E° ) (E° - E° ) (E° - E° m ) 
u v n n' 7 n n" 7 n n' u/ 


+ 


5 h5!1 


5 H. ( P. 


24 ( E n- E n.) K " E n"> <<»■<«.) ^ < E X«> (E n’- E n”> 


+ 


1 H 


(!) 


1 h/2 


6 <«.») (E n--<".> K~ E °n^ 6 (E n- E n'” } 


1 h/!> 


i h/. 1 ) 


6 (<-<.) <«-<,„) 6 <<„-<•»> 


1 


1 H<H 


24 <*£" V ^ E n'- E n") < E X«.> <<“<■> (E n'- E n"> 


1 H 


(*) 


1 H. ( ?> 


24 n 0 On 0° 24 0 0 00 , 0 0 

(E^E n ,„) (E n „-E“,„) (<-V (E n -E n ,„, (E n „-E n ,„) (E^-E^J 


1 


i h?;>. 


M (<-<■"> (<-<■> 24 (<-<"•> (E n-<'> ( <'- E n"' 


i «<:>. 


i h/JI 


24 0 0 0 0 0 

^ (E°-E t )(E -E mt ) (E -E n ) 

n n T/ n n MT/ ' n f n n 


24 (E° -E 0 ,) (E 0 -E° im ) (E u -E u ) 
n n' n n ,M n n n'" 


r(l) 


* This form will be denoted by H ... in this part of the discussion to make 
things easier to write. 
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i h. ( !1 


1 


24 00 00 00 24 00 00 00 

(E -E „)(E -E „)(E -E ,„) (E -E ) (E ,-E „) (E ,,-E ,„) 

' n n" n n* n"' n n m n n n/ n’ n" n" n" 1 


1 

24 o o oo o o 

(E -E )(E -E WE -E ) 

1 n n nM n» n" ?M n n T; 


JL H??l 

24 0 0 00 00 

(E -E ) (E -E ) (E -E ) 
1 n n" M n" n ,,,M n» 


Some of the above tactors required two index changes to bring them into the 
desired order. For example, in the last factor above, n ,M and n" were inter- 
changed first, then n" and n'. To combine the above factors, it will be neces- 
sary to identify each factor by "factor number" 1, 2, ... depending on the order 
in which they are found. With this in mind, one can recombine the above factors, 
as follows: 

Factor numbers (4) and (7) give 

1 1 _ 1 \ _ -H.ty 

6 (E° -E° ) (E° -E° ) E° -E* E°-E ° J 6 (E° tI -E° MI ) (E° -E° „,) (E° -E° ,) 

n' n n,/ n" n"' n n"' n n" n" n"' n n'" n n' 


This result then combines with factor numbers ( 10) and (12) , giving 


24 o o 0 0 00 

(E -E° ) (E -E„ t )(E ,,-E „ ) 
n n' n n"’ n" n"' 


Factor numbers (5) and (6) give 


1 H.ty / 1 _ 1 

6 (E 0 ,-E°,„) (E° r E°„) \E°-E°, E^E 0 ,,, 
n T n ,M/ n' n" ' n n' n n 1 " 


00 0 „o 00 

6(E ,-E ) (E -E ) (E -E ,„) 

n' n" 7 n n' n n ,M/ 


This result, in turn, combines with factor numbers (13) and (14) to give 


24 (E°-E° ) (E°-E 0 ) (E° -E° ) 

n n* n n M,/ n’ n" 


49 



Factor numbers ( 8) and ( 16) give 


24 


h5?1 


IE°-E°, 
n n' 


o o 
E - E 
n n'" 


H (i) 

Xl« , , 

24(E°-E° fII ) (E°-E° ) (E° -E° ,) 
v n n ,M n n" ' n n 1 

(C-) 


Factor numbers (15) and (18) combine to give 

1 1 


1 h/?. 


h/?. 


24 


<<-*;.) \ E n- E n>. (eO-E^^eO-E^^eO-eO,,) 


Factor numbers (11) and (19) give 

1 


1 H.^il 


24 


E°-E° 


n n 


in 


E°-E°, 
n n' 


(D') 




n 0 0 o 0 0 

24 (E° -E „)(E-E ) (E -E,,,) 

n' n" n n" 7 n n m 


This result, in turn, combines with factor number (3) to give 


6 h5!1 


24 (E° - E° ) (E°-E° ) (E° -E° ) 

n n" n n'" n* n" 


(E») 


This result then combines with (B') to give 


6 H. ( !l 


24 (E°-E° ) (E°-E° ) (E° - E° ) 
' n n' 7 n n n/ n n m 


(F-) 


Factor numbers (9) and (17) combine to give 

i h. ( ! ) . 


i h.<!>. 


24 0 0 OOlOn 00l?4 0 0 00 n0. 

(E ,,-E ,„)(E ,-E „) \E -E , E -E „/ Z4 (E ,,-E ,,.) (E -E ,) (E°-E „ 

' n" n"’ 7 n’ n" 7 n n 1 n n" n" n'" n n' 7 n n" 


This result combines with factor number (2) to give 
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6 h5!1 


24 0 0 0 o o o 

(E -E „,)(E -E ) (E -E „) 
' n" n"’ 7 a n’ 7 n n" 


(G 1 ) 


Combining (G’) with (A’) one obtains 


6 


24 (E 0 -E 0 | )(E 0 -E 0 , i )(E 0 -E°, m ) 
' n n’ 7 n n" n n"’ 


(H') 


Finally, combining factor number (i) with (C 1 ), (D’), (F’) , and (H’) one 
obtains 


10 1 1 6 6 y 1 V 

24 + 24 + 24 + 24 + 24* h .. h 


2 


H (‘) H <’> H (1> H (1) 

nn f n f n" n"n m n™n 


o o 


.o o 


o 0 


n ! ^n n n ^n,n f n Mt ^n,n T ,n !T (E -E ) (E -E ) (E -E ) 
* ’ n n !/ n n M/ ' n n m/ 


Note that all nineteen terms are accounted for and that the five triple sums in E 
have now been replaced by 


( 4 > 


n 


2 2 


2 


,(•) „(>) „(<) „(*) 

nn’ n’n" n M n"* n"’n 


n’^n n’ ,: *n,n’ n"' : J t n,n',n" (E°-E ) (E°-E° ) (E°-E° ) 

» ’ ' n n' n n" 7 n n ,M7 


(5T) 


plus terms (ACE) and (BD). 

( 4 ) 

An additional simplification in E^ is possible by grouping the three 

double sums containing combinations of H ^ and H ^ matrix elements. To 

n’n n’n 

avoid having to rewrite the terms in question, they will be referred to as 
"first", "second", and "third" terms depending on the order in which they 

appear in E^ . Adding the first part of the first term to the first part of 
the second term results in 


1 

3 


2 

n’^n 


2 


nn’ 


n"^n,n’ 


h *!\, h *h + H<1) „ 

n’n" n"n nn" n"n* n’n 

(E° - E 0 ,) 
n n 7 



Then, by interchanging n" and n* in one of the factors, one gets 
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I z z 

n'^ n"^n, n' 


( 1 ) ( 1 ) ( 2 ) 
nn' n'n" n"n 


E° -E° 
n n* 


E°-E° 
n n" 


E u -E 


n' n" / 


1 


E°-E° 
n n" 


/ 

1 

1 n 

E°-E°. 

E° -E° ( 

\ n 

n 1 

P 

P 


. H (1) , H^,, H (2) f(E° + E°, - 2E 0 ,,) + (2E°, - E°. - E u )] 

_i ^ nn 1 n'n" n"n n n' n" ; v n' n" n 7 ' 

3 n'^n n”^n,n' (E°-E° t ) (E°-E°„)\ (E° - E° ) 

* v n n ,M n n" 7 ^ K n » n "' 


H (1) , H ( J } n H (2) 
_ ^ 1111 n n n”n 


n'^n n”*n,n' (E° - E°.) (E° - E°„) 
n n' 7 n n” 7 


Next, add the second factor of the first term to the second factor of the second 
term in the same way as above. This sum is 


h (2) , h^,, h<;> 

nn T n T n n n n n 


y y 

l— I £—J On 0 0 

n'*n n”#n,n' (E - E°,)(E - E ) 
n n T,v n n" 7 


By changing n" to n' in the second factor of the third term, one obtains a 
similar expression. Combining these three terms, one obtains 


Z Z 

n'^n n”^n, n' 


(H (2) , H<J> n 

' nn' n'n" 


H 


<*) +H (l) „(*) „(D 


n"n 


nn' 


H' 


n'n' 


H 


+ 

n"n nn' n'n” 


H (2) ) 
n"n 7 


(E° - E° ) ( E° - E° ) 
n n' 7 n n" 7 


(3T) 


(4) 

The last set of terms that need to be simplified in E^ are those con- 
taining products of the form H^, H^„ . The method is the same as 

& nn' n'n” n"n 

above: One takes the terms in the order in which they appear and interchanges 
n" and n* in the second factors of each term. This gives the following expression: 
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(!) 


M) _„<*)• 


, , , ,v r2(H'/ -H w ) 2(H ,,,,-H ) 

i y V H (l) H (i) H (l) I n n * nn ; f n"n" nn 

3 ,V ..jA 7 . nn' n’n" n"n \ ,_o _.o .2,^1° t^o * /T .o t-,0 , 2,t-,° t .0 

n'^n n"^n,n' ] (E -E ,)4E -E „) (E -E )^(E -E , 

L n n' 7 n n" 7 n n" 7 n n' 


(H ( J\ - H (l) ) 
n f n T nn 7 


(H^ „ - H (1) ) 
n"n n nn 7 


OOoO 0 OOoOO 

(E -E ,) 2 (E ,-E „) (E -E ) 2 ( E ,-E° ) 
n n* 7 n* n TT/ « n H/ r» T n |T/ 


0 

5 -E 
n n' 


n' n' 


(H ( ^ - H ( !\) 

n"n" n'n' 


(H ( ^ - H ( I\) 

n"n" n'n' 


0 0 0 0 o 

(E -E ,)(E° -E .. ) 2 
n n' 7V n' n" 7 


(E°-E°„) (E° ,-E° ) 2 
n n" 7 n* n" 7 


Rearranging the brace, one separates the terms as follows 


H 


A) 

nn 


0 0 9 

(E -E y 
n n 7 


2 1 
4 - 




e°-e° m e^-e 0 , 

V n n" n' n' 


(E ~E „)' 
n n" 7 


E -E , 
n n' 


E n' E n" / 


+ H 


(i) 

n'n' 


i 

2_ + i ) 

1 

1 1 ) 


(E° - E 0 ,) 2 
n n' 

1 E°-E° E° -E° i 

1 n n" n' n"J 

+ (E n'- E n") 2 

1 E 0 -E 0 „ e -e 1 
\ n n" n n ') 



+ H 


( 1 ) 

n"n' 


(e°-e 0 m ) : 

n n" 7 


F°-F° E° - E° 
n n' 11 n 


(E° ,-E° ) 2 1 E°-E° E 0 -E° I( 

n' n" \ n n' n n" 


- H 


(4 

nn 


|(E°-E° t )( E 0 -E° t ) (E° -E ,,) 
In n TI v ix n ?; n ? n" 


2(E° ,-E° )+( E° -E° ) 2(E 0 ,-E 0 |t )-(E 0 -E 0 T )1 

n' n" n n" 7 n' n" n n T 

.0 -.0 > , T ^0 ^0 


(E - E ) 
n n' 7 


(E - E ) 
n n" 7 


+ H<| ) | [ 2(E <> ,-E° ) + (E°-E° ) - (E°-E° )] 
n'n' n' n” 7 n n" 7 n n ,7J 


(E°-E°,) 2 (E° -E° ) (E° -E° ) 
n n' n n" 7 n' n" 7 
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+ H 


(i) 


"n» 


n"n' 


[2(E»,- E »„) ■<■<,) + (<-<■, » ~j 

(Ej-E*,) (E»-Ep 2 (E«,-E»„) / 


With a little extra effort on the first parenthesis, one simplifies the brace as 
follows: 



3 H (1) (2E°-E° -E° ) 
nn n n 1 n" 7 

(E°-E° ) 2 (E°-E°,) 2 

v n n* !/ ' n tv 7 


Thus, the group of terms just considered equals 


E E 


H (, ». H<‘> 


H 


(i) 


nn’ n'n" n"n 


H 


(i) 


H 


<*) 


n'*n n"*n,n» (E°-E 0 .) (E 0 -E°„) 
n n' n n" 


oo oo 
E -E „ E -E , 
n n" n n' 


- H 


(i) 

nn 


E E 

n'#n n"#n, n' 


(2E° - E° - E° ) 
nn' n'n" n"n 1 n n' _n" ; 

(E° - E°,) 2 (E° - E° ) 2 
' n n' 7 n n" 


(4T) 


( 4 ) 

At this point, it is best to write down E fl incorporating all the simplifi- 
cations just obtained and the additional terms generated. One obtains for E^ 
the following (continued on next page) : n 


54 



f H ( III 2 ) 

e (4 > = h ( 4 > + y < rm ' ^ 

n 1111 n'*n \(E° - E ,) 
v ' n n' 




nn' n'n mr n'n 

(E° - E°,) 
n n' 7 


H (l) , (H ( J\ - H (1) ) 2 

nn 1 n'n n'n* nn 

(E° - E°,) 3 
' n n' 7 


+ H (2) , H^) (H<I>. - H (1) 
' nn 1 n'n nn 1 n'n n'n 1 nn 

(E° - E°,) 2 
' n n* 7 


H (1) , H ( J ) (E^-E^) (H ( J\-H (2) ) 

nn' n'n n' n 7 nn' n'n n'n* nn 
■t* "f* 


2 (E° - E° t ) 2 
' n n' 7 


2 (E° - E° ) 2 
' n n' 7 


Z Z 


„(i) „<■> H (‘) 

nn' n'n" n"n 


n'#n n"*n,n' (E°-E°,) (E° - E° n ) 
’ ' n n' 7 ' n n" 7 


H 


( 1 ) 


"n" 


n"n' 


H 


(0 

n'n' 


\ E° - E° 
n n' 1 


E° - E°, 
n n' 


- H 


(i) 

nn 


H (1) , H^,, (2E° - E° - E° ) 

nn' n'n" n"n n n' n" 7 

0 0 O ft 0 9 

< E n- E n'> < E °- E »") 


+ z z 

n'^n n'^n, n' 


Z Z 

n'^n n"#n,n' 

(h (2) , h ( J\, h (4) +h (1) , h (2) „ h ( 4) +H (I) , H<J\, H<«) 
' nn 1 n T n !T n M n nn* n’n" n t! n nn T n T n n n"n 


z z z 


( E° - E 0 ,) (E° - E° ) 
' n n' 7 ' n n" 7 

(i) (i) (l) (l) 

nn' n'n" n"n"' n"'n 


n'*n n"*n,n* n"'#n,n',n" (E°-E° t ) (E°-E°„) (E 0 -E°„ t ) 

' n n f/ n n ,T/ ' n n Tt,/ 


A z z 


to (i) (j> (i) 

nn' n'n" n"n' n'n 


I i 


2 LJ LJ 0 0 0 0 

n'^n n'^n, n* (E -E ,) (E -E ,,) 
’ v n n' 7 ' n n" 7 


E° -E° 

\ n n' 


E -E „ 
n n" / 


-7 Z Z 


in (i) («) „(D 

rl tt xl xl , xl . 

nn M n T! n nn f n’n 


4 LJ LJ 0 0 0 0 

n'* n n”*n,n' (E -E ,) (E-E ,,) 
x n n f 1 n n/ 


n n T 


E°-E° t 
n n* 


E° -E° 
n n" 
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The last two terms correspond to terms (ACE) and (BD). These terms will 
combine with others generated by increasing the range on n" and n m later on. 

Since the result for E is now expressed in terms of denominators 

n 

which exclude factors of the form (E^-E^,,) • (E®,-E^ m ) * anc * ( E n" _E n'"^ * 

one may now put back those terms which are not affected by this type of 
singularity. The last two terms do not fit this description, but as was pointed 
out earlier, they will combine properly when all the modifications are included. 

One may now express in terms of sums having the following 

form:* 

2 « • 2 2 q . 2 2 2 q • 

n'An n’^n n"=£n n'^n n"An n"'An 


Thus to change the double sums one needs to take 


2 2 <5= 2 2 <3- 2 Q(n" = n'>- 

n'^n n"An,n' n'An n'^n n'^n 

( 4 ) 

With this in mind, one replaces the triple sum term in E^ by 


(l) (i) (1) (i) 

XT , fl , XX t|f XX ... 

nn T n T n n n n n™ n m n 


Z/ Z1 Z/ on 00„0„0 

n T ^n n M ^n n MT: £n (E -E° ) ( E -E )(E -E ) 
v n n T/ n n M/ n n ,M/ 


- 2 2 


(H^ 1 * H^„ , H^) 

v nn ! n T n n n"n T n'n' 


n’^n n"An, n' (E°-E^,) 2 ( E °- E !„) 


n ir 


nn n f n M n"n 


n"n" 


H 


( 1 ) 


n*n 


- 7 7 

n'^n n"^n,n' (E°-E° ) (E 0 -E°„) \e° -E°, t E°-E°, / 

v n n' 7 n n"' n n" n n' 


.o „o 


0 ^0 


(v) 


* Let Q denote any given expression. 
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The newly generated double sums are incorporated into the existing terms in 

before the above modification of the double sums is performed. In doing 

n (4) 

this one notes that the last term in (v) cancels one of the terms in E , the 

(4) n 

second term in (v) adds to the next to the last term in E . This last term will 

n 

(4) 

be incorporated later on. In referring to , one notes that the next term to be 
simplified is the second double sum. When this term is modified, one obtains 

H (1) , H^,, (2E° -E° -E° ) 

- H ^ V y n n n'n" n"n v n n' n n/ 

1111 n'^n n"#n (E°-E 0 ,) 2 (E°-E° n ) 2 

n n' 7 ' n n" 7 


+ 2 H' 


(H (l) , H^) H^, 

Z nn’ n'n 7 n’n* 

o o <t 

n'*n (E - E ,) 3 

' n -n i ' 


The new factor generated above cancels the middle factor obtained when the 

fourth term in E^ is multiplied out.* The fourth term (a double sum) of E^ 
n n 

that requires modification becomes 


Z Z 

n'^n n"#n 


(H (2) , H ( J\, + H (1) , H ( f H ( ^ + H (1) t H (2) ) 

' nn' n'n'-' n"n nn' n’n” n"n nn' n'n" n'T 

(E° - E° f ) ( E° - E° ) 
n n' 7 n n" 


Z 

n'^n 


+ H (2) H^) 

nn' n'n nn 1 n'n 7 n’n* 

(E° - E 0 ,) 2 
n n' 7 


(H (1) , H^) H (2) , 
y nn' n'n 7 n'n* 

Lt o o , 

n'^n (E - E ,) 2 

n n' 7 


Here again, one notes that the new terms generated above cancel some of the 
(4) 

existing terms in E^ . 

Before incorporating all the terms generated one should take care of the 
last two terms in E^ . Combining the first term with the middle factor of (v), 
one obtains 


The last form of E ' 7 is referred to in each of the cases that follow, 
n 
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-5 2 2 


H (1) . H (1) 


nn' 


H 


(i) H (i) 


n'n" n"n' n'n 


2 t-i 00,0 o 

n'^n n'^^n' (E - E ,) 2 (E , - E ,,) 
n n' 7 ' n* n" 7 


Now, if one puts the factor corresponding to n M = n, back into the above expres 
sion, one obtains 


-I 2 2 


H (l) H (l) H (l) H (l) 


2 n'*n n"£n’ (E°-E°,) 2 (E 0 | -E°„) 
n n 17 n’ n" 


'n 1^ y 

” 2 h 


(H (1) , H ( |V 
rm T n , n / 


2 n'^n (E°-E° ) 3 
n n' 7 


, (y) 


(both signs are negative because, when n" = n, (E^, - E^ n ) is negative). When 
the first term in (y) is rearranged, the result is 

,/ (1 > 

V - 

2 


2 


nir 


2 


H (1) H 10 

n'n" n"n' 


H 


( 1 ) 

nn' 


o o l *_/ oo 

n'*n (E -E ,) \n"*n' (E - E „) 

' n n' 7 n' n" 7 


(E°-E° ) 
' n n* 7 


Now, by using equation (26) in the form 


(E^ - H ( *\) = V 
' n 1 n’n’ 7 .H 


m (1) H (l) > 
1 n’n" n"n 17 


0 0 

n^n’ (E f - E „) 

n T tv ' 7 


n r 


(y) becomes 


-f 2 


H (1) , (E^ - E (2) ) 

nn 1 n T n n 1 n 7 


n 


n'^n 


o 0 O 

(E - E ,) 2 
n n' 7 


-I 2 


(H (1) , H ( J } ) 2 
nn* n’n 


2 0 0 3 
n'#n (E -E 0 ,) 3 
n n' 7 


(y T ) 


By performing the above manipulations, one eliminates what appeared to be a 

(4) 

"bad" term. The same analysis takes care of the last term in E^ . Taking out 
the factor corresponding to n" * n', one obtains, for this last factor, the expression 
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Now, by separating out terms in the first factor above and by using equation 
(26) , one obtains 


7 Z - " 


(E (2) - H (2) ) H (l) , 

' n nn 7 nn' n’n 1 


n’^n 


(E° - E 0 ,) 2 
n n' 7 


n"#n 


H (1) „ (E (2) -H (2) ) 

nn” n"n v n nn 

(E° - E°„) 2 
n n" 


Let n" = n' in the second term above, and equation (z) becomes 


-i Z 


H (1) , (E (2) - H (2) ) 

nn 1 n'n x n nn 7 

(E° - E 0 ,) 2 
n n' 7 


(H (1) , H ( JV 

Z nn 1 n'n 7 

o n i 
n'#n2(E -E u ,) 3 
' n n 1 7 


By combining (y 1 ) and (z 1 ) and all the other terms generated, one finally obtains 
(continued on next page) : 
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e ( 4 ) =h (4) + y - — 11 11 ^ 

n nn o 

n'^n (E - E .) 

n n' 7 


/u(*) tr(3) + ^ ) 

(H nn' H n'n H nn' H n 'n J 


h (2) , H ( ^ 

V nn 1 n'n 

Li oo 

n'^n (E -E ) 
n n' 7 


H (l) , H (l) , [ (H ^,) 2 + (H (1) ) 2 ] 

e ( 2 ) y mi 1 n'n + y nn 1 n'n 1 n' n' 7 nn 7 1 

n n'*n (E° - E°,) 2 n'^n (E° - E°,) 3 

' n n' 7 n n' 7 


(H^ H*J\, H^ 4) ) (2E° -E° -E° ) 
jj(i) V y nn* n'n" n"n ; n n' n"' 

nn n'*n n"*n (E° - E 0 ,) 2 (E° - E 0 ,,) 2 

n n' 7 n n" 7 


E E 

n'^n n'^n 


(H ( 2 ) ,H ( J ) „H ( ‘ ) 

v nn' n'n" n"n nn' n'n" n"n nn' n'n" n"n 


(E° - E°,) ( E° - E° ) 
n n' n n" 


(H (1) , H^,, H (4) H<J> ) 

Y> y y nn' n'n" n"n'" n"'n 

n'^n n"*n n"**n (E° - E° f ) (E° - E°„) (E° - E°„,) 

n n r n n" 7 n n'" 


This answer may now be used for numerical calculations requiring perturba- 
tion terms up to fourth order. With the aid of computers this type of calcula- 
tion is becoming more prevalent. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
November 1965 
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